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Abstract 

For fields that vary slowly on the scale of the lightest mass the logarithm of the 
vacuum functional can be expanded as a sum of local functionals. For Yang-Mills 
theory the leading term in the expansion dominates large distance effects and leads 
to an area law for the Wilson loop. However, this expansion cannot be expected 
to converge for fields that vary more rapidly. By studying the analyticity of the 
vacuum functional under scale transformations we show how to re-sum this series 
so as to reconstruct the vacuum functional for arbitrary fields. 



1 Introduction 



When the Schrodinger representation vacuum functional of a massive quantum field theory 
is evaluated for slowly varying fields its logarithm reduces to a sum of local integrals. For 
example, the vacuum functional of a free scalar field theory in D + 1 dimensions with 
mass m is ^[<p] = {<p\0) = exp —\j <i D x (fy/— V 2 + m? ip, so that if the Fourier transform 
of ip vanishes for momenta with magnitude greater than the mass, m, the logarithm of \I/ 
can be expanded in the convergent series 

The terms of this expansion are local in the sense that they involve the field and a finite 
number of its derivatives at the same spatial point. The same is true for an interacting 
theory in which the lightest particle has non-zero mass, because the logarithm of the 
vacuum functional is a sum of connected Feynman diagrams, and since massive propaga- 
tors are exponentially damped at large distances these diagrams reduce to local integrals 
for slowly varying fields. In the case of 3+1 dimensional Yang-Mills theory the lightest 
glueball mass must depend non-perturbatively on the coupling so that the local expan- 
sion of the vacuum functional could not emerge from the usual semi-classical expansion 
of the functional integral. Nonetheless we would expect to obtain for slowly varying 
gauge-potential 

W[A}=ln^[A] = [(FxfaitrB-B/m + aQtrD AB-D AB/m 3 + 

a 3 trB- (B A B)/m 3 + a 4 tr B ■ B B ■ B/m 5 + ..) (2) 

where B is the Yang-Mills magnetic field V A A + A A A, and D is the gauge covariant 
derivative. (This is qualitatitively different to the case of pure Abelian gauge theory where 
W = — / <^ 3 x d 3 y B(x) -B(y) / (x— y) 2 is conformally invariant and cannot be expanded 
as an integral of local quantities.) The unknown coefficents a, are dimensionless constants. 
In principle these coefficents can be determined from the Schrodinger equation. In Jl| we 
showed that the Schrodinger equation does not take its usual form because the process 
of removing the cut-off does not commute with expanding in terms of local quantities. 
However, by re-summing in terms of the cut-off an appropriate Schrodinger equation can 
be constructed, and we suggested a method of solution that does not rely on a semi- 
classical perturbation theory. Alternatively the first few terms can be computed in lattice 
gauge therory. Approximating the vacuum functional by the first term of (0) was originally 
suggested by Greensite @, who used Monte Carlo techniques to compute its coefficent 
for the gauge group SU(2) in 2+1 dimensions |J and, with Iwasaki, in 3+1 dimensions 
||]. When these results are translated from lattice to continuum fields they give ai/m — 
fexp [^-7r 2 (/? — Po)], [|J, where a is the lattice spacing, (3 the conventional coupling that 
enters the Wilson action for the 3+1 dimensional theory, and (3 « 1.74. The 2+1 
dimensional result was extended by Arisue in |J, and the case of gauge group SU(3) 
considered by Chen et al in J/J. For SU(2) the 2+1 dimensional theory yields W^A] = 
Jd 2 x (!±trB 2 + ^trD i BD i B + ..y §, where b x = 0.91 + 0.02, b 2 = -0.19 + 0.05 and e 
is the continuum coupling constant in 2+1 dimensions corresponding to the normalisation 
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of the Lie algebra generators such that trT A T B = —5 AB /2. (Note, that whilst the naive 
continuum limit of the lattice strong coupling expansion || would yield a similar expansion 
the coefh cents would be wrong as they do not scale appropriately in the continuum limit.) 
To get some idea of the relevance of these results to large distance effects we can attempt 
to compute the string tension. We will simplify the calculation by retaining only the 
leading term. In the 2+1 dimensional theory this leads to an area law for the Wilson loop 
H because the leading term in W3 is the action for two-dimensional Yang-Mills which is 
free and confining. Applying the calculation given in || to the result of || gives the string 
tension as 03 = 3e 4 /(16&i) ~ 0.21 e 4 which should be compared with the recent direct 



Monte Carlo estimate a = (0.112 ± 0.001) e 4 [10]. Similarly the large Wilson loops in the 



3+1 dimensional theory may be estimated using just the first term in (|]) which reduces 
to a calculation in the three-dimensional theory, and thence to a calculation in the two- 
dimensional theory, so that the area law emerges from a kind of dimensional reduction 



llj. Putting this together gives the string tension in the 3+1 dimensional theory as 
a 4 « 0.21e 4 with e 2 = m/{2a 1 ) so that a 4 a 2 « 3.4 exp [-^((3 - 1.74)]. For (3 = 2.85 this 
yields a 4 a 2 « 0.0086 which compares with the accurate direct Monte Carlo measurement 
a 4 a 2 = 0.00363(17) HT2] . Some comments are in order. Firstly, the coefficents a±, b\, b 2 are 
obtained by applying the naive continuum limit to a fit of the lattice vacuum functional 
obtained for slowly varying plaquette variables. This is justified since the transition 
from plaquette variables to continuum ones is itself a derivative expansion in continuum 
variables. If we had computed Wilson loops using the form of the vacuum functional 
originally found in terms of plaquette variables we would also have obtained an area law, 
but with a different string tension that would not scale appropriately in the continuum 
limit. This is because the integration over plaquette variables would include rapidly 
varying fields that were excluded when we made the transition to continuum variables, 
and for which the form of the vacuum functional used is inappropriate. The calculation of 
the string tension using just ai and bi, although numerically incorrect by about a factor 
of two, does at least scale properly. Presumably the numerical value would be improved 
were we to include the other terms in the expansion of W. In spite of this qualified success 
in describing large Wilson loops we can only expect the local expansion to converge for 
configurations that vary slowly on the scale of the mass of the lightest glueball so it would 
appear not to be relevant, for example, to the computation of the glueball spectrum as this 
involves heavier particles. The purpose of this letter is to show that in fact the vacuum 
functional for arbitrary A can be reconstructed from this local expansion (Q), given the 
coefficents cij. 

We will work with the functional integral representation 

tf[A] = JvAe- s[A] - s * [A ' A] (3) 

where S[A] is the Yang-Mills action gauge-fixed in the Weyl-gauge A = 0. Space-time is 
Euclidean with co-ordinates (x, t) and t < 0, so 

S[A] = -jj J d 3 x dt tr [A 2 + (V A A + A A A) 2 ) (4) 
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The boundary term in the action is 



S b [A,A] = -j? Jd^tr [{A-A)-A)\ t=0 (5) 

The boundary value of A is to be freely integrated over, i.e. we will not impose a condition 
such as A(x, 0) = 0. We will assume that at spatial infinity the source A is a pure gauge 
A ~ g(x.)~ 1 Vg('x) . St, is chosen so that \I/[A] is invariant under the gauge transformation 
5^ A = Vuj + [A, u], since the effect of varying the source A may be compensated by 
gauge transforming A. As u cannot depend on time this is the residual gauge symmetry 
of S[A] that preserves the gauge condition. Functionally differentiating with respect to the 
source leads to an insertion of A. Given the Wick rotation between Euclidean time, t, and 
Minkowskian time, this yields the Schrodinger representation of the canonical momentum, 
which in this case is the non-Abelian electric field, E = —ig 2 5/5A. 

Consider the effect of a scale transformation on the configuration A given by 

AS W = ^A(^x). ( 6 ) 

For small s the scaled field is trivial everywhere except the vicinity of the origin, since 
for x = we have A s (0) = s~ 1/2 A(0), but for x ^ A s (x) « s -1/2 -1 (x)V0(x). For 
large s the field varies slowly in space since now A s (x) m s _1 / 2 A(0). By studying the 
analyticity of ^[A 3 ] in s we will be able to express its value at s — 1 in terms of its value 
for large s, where we can apply the local expansion (0), and its value for small s, which 
is reliably computed in the usual semi-classical perturbation theory and turns out to be 
negligible. More specifically, we will show that \1/[A S ] can be analytically continued to 
the complex plane with the negative real axis removed. This enables us to compute the 
contour integral 

in two ways. We take C to be a key-hole shaped contour running just under the negative 
real axis up to s = 1 — s Q , around the circle of radius s centred on s = 1 and then back 
to s — — oo running just above the negative real axis. If we take sq to be large then we 
can compute the integral using the local expansion for W. Each term in this expansion 
scales so we can express it in terms of A rather than A s , 

W[A S ] = J d 3 x(a lS - 1/2 trB-B/m + a 2 s- 3/2 trD AB-D AB/m 3 + 

a 3 s~ 3/2 trB- (BAB)/m 3 + a 4 s" 5/2 trB-BB-B/m 5 + ..) (8) 

This yields an expansion of \I/[A S ] in inverse powers of s — 1, with coefncents that depend 
on the original configuration, \1/[A S ] ~ I] (s — l) _n ^ n [A], enabling us to compute /(A) as 

We can also evaluate the integral by collapsing the contour C until it breaks into two 
disconnected pieces, a small circle centred on s = 1 and a contour that just surrounds 
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the negative real axis. The integral over the circle gives ^[A]. By taking A to be real, 
positive and very large the contribution from the negative real axis will be exponentially 
suppressed, (provided it is not singular, which we check in perturbation theory), so we 
obtain for large A 

*[A1 ^'> [A | (10) 

which provides a re-summation of the local expansion. Note that only terms of order 
up to s~ n in (H) will contribute to ip n . We might expect to obtain an approximation by 
truncating ( |10|) at some order in A. 

As an illustration we will show how the vacuum functional of a free massive field 
theory may be reconstructed from its local expansion. First scale the field by setting 
(p s (x) = s~ ( - D ~ 1 ^ 2 ip(x/y/s). The vacuum functional for the scaled field is then exp W[p s ] = 
exp — \ j d D x p\/ — V 2 + sm 2 <p which can be continued to an analytic function in the 
complex s-plane with the negative real axis removed, and is finite at the origin. For large 
s we can expand this, or alternatively ([!]) evaluated for ip s , in inverse powers of s — 1 
obtaining the local series 

which yields an expansion of $ in inverse powers of s — 1. The term quadratic in (p is 
just ly^ 5 ] itself, so if we concentrate on this term rather than the whole of \l/ we should 
consider 

2mJcs-l [y 1 iy/n^nl (n-l/2)7 ^\ m 2 J y v ; 

The integrals will exist for all n provided tp has a momentum cut-off, ko say, but the 
integral will converge for all A and ko because of the n\ in the denominator. This is in 
contrast to our original expansion (0) which only converges for ko < m. The series (|T^ ) 
results from expanding in powers of A the exponentials in 



"' r d D Xip | —— e - A ( 1 " v2 / m2 ) _|_ r dX J- e - x ^ 2 ^ 2 Ul-^-)) V , (13) 



A 




2V? J Y Jo y/\ 

which is just 

-U^V=v^ v + ^l^ (J" d x-L e -*-^) „ (u) 

We see that as A — > oo the series tends to VT[<^], as it should. Furthermore, if we keep A 
large, but finite, the error in approximating W[ip] by the series (|T^), as given by the last 
integral in (|T4]), is exponentially suppressed. For a finite value of A we can also truncate 
the alternating series (O) at order A™ with an error smaller than the first neglected term, 



so that by taking n sufficently large in comparison to A and k this error can be made 
small. In conclusion, we can represent W[ip] for cut-off, but otherwise arbitrary, ip, by the 
local series flT2D truncated at order A n where A and n are both large and chosen to give 
acceptable error. 
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2 Analyticity of the Vacuum Functional 



To construct an analytic continuation of \1/[A S ] we adopt a similar approach to that used 
in complicated by having to work with three spatial dimensions. We consider scaling 
each dimension separately, so we set A s = (sis 2 S3)~ 1 ^ 6 A(x 1 / \/s v x 2 / x 3 / y/s 3 ) and 
show that \I/[A S ] is analytic in si,s 2 ,s% separately. Firstly we interchange the names of 
the Euclidean time, t, and one of the spatial co-ordinates, x 1 say, in the functional integral 
(|3|) which we now interpret as the Euclidean time-ordered vacuum expectation value 

^[A s ] = T(0 r | exp (jj J dx 2 dx 3 dttr (A s • A')^j ^ |0 r ) (15) 

where |0 r ) is the vacuum for the Yang-Mills Hamiltonian defined on the space x 1 < 
in the axial gauge A\ — with a boundary term in the action \ J dx 2 dx 3 dttr (A ■ A'), 
where the / denotes differentiation with respect to x 1 . Expanding the exponential gives 

*[A S ] = 

f dt n [" 1 dtn-i.. / 3 dt 2 f 2 dti J] f— 3 / dx 2 dx 3 A s Ri (U, x 2 } x 3 )) 

n J —oo J— oo J— oo J— oo \ y J J 

(0 r \A Rn (0, x 2 n , x 3 n ) e^-^ H ..A' R2 (0, 4 4) e {t2 ~ ti)H A' Rl (0, x 3 ) |0 r ). (16) 

Here Ri stands for both Lie algebra and spatial indices. The time integrals may be done 
after Fourier transforming the sources. To do this we define the Si-independent Fourier 
mode 

a(k,x 2 ,x 3 ) = (s 2 s 3 )- 1 ^J dte- ikt A(t,x 2 /^ 2 ,x 3 /^ 3 ), (17) 

so that 

A s {t,x 2 ,x 3 ) = ±J dke lkt/ ^ s s- 1/6 a{k,x 2 ,x 3 ). (18) 
Substituting this into ([IB]) gives 
#[A S ] = 

Hft [ dhdx 2 dx 3 s\ /3 a Ri (ki,x 2 ,x 3 )) 6 

n i=i v J / 

{0AA '^°' x « * AH-'ig-'V ^^ (0 ' I? ' I?) |0r> ' (19) 

This makes explicit the Si-dependence of \I/[A S ]. Although s\ was originally real and 
positive we may use this expression to define an analytic continuation to complex values, 
yielding a function that is analytic away from the zeroes of the denominators in (|T9|). 
Since the eigenvalues of the Hermitian Hamiltonian are real, the singularities lie on the 
negative real Si-axis. Similarly we may show that ^[A 55 ] continues to an analytic function 
in s 2 and s 3 and, by setting S\ = s 2 = s 3 = s, that ^[A 5 ] continues to an analytic function 
in s on the complex plane with the negative real axis removed. 

We have seen that for small values of s the configuration A s is non-trivial only over 
a short distance about the origin. Since asymptotic freedom under-writes semi-classical 
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perturbation theory at short distances we can reliably use this to calculate \1/[A S ] for 
small s. We will assume that Symanzik's work on the vacuum functional of </> 4 theory ]nj 
can be generalised to the present case, and further, that gauge invariance implies that 
the source A needs no renormalisation WM. Thus we assume that the cut-off necessary 



to define (|3|) may be removed leaving a finite vacuum functional that depends on the 
source A, an arbitrary mass-scale \i and a renormalised coupling, g[p). (This assumption 
is supported by the scaleing behaviour observed in lattice estimates of \V, ||, M, ||). 
Now if in computing \I/[A S ] we were to choose a new unit of length so as to undo the 
scaleing of A, and at the same time we scale \x appropriately, then nothing would change 
so we must have that \I/[A S ] computed using /i and g(p) is equal to \I/[A, ] computed using 
fJi/y/s and g(n/y/s). As s decreases so does the coupling, since in perturbation theory 
g~ z / ' y/s) = g~ 2 (fi>) — UN (Ins) /(A8tt) for gauge group SU(N), so for small enough s 
it is sufficent to take the tree-level approximation 

MA] PS e -S[A cl ]-S b [A cl ,A] ^20) 

where A c i satisfies the Euler-Lagrange equation 5(5 , [^4] — Sb[A, A]) = under variations 
of A that are arbitrary at the boundary t — 0. Now 

6 U J d 3 xdttr [A 2 + (V A A + A A A) 2 ) + j d?*tr ((A - A) ■ A) \ t=0 ) = 

J d 3 xdttr ((AdA) -5A[a + VAB + AAB + B a a)) 

+ J rf 3 xtr ((A- A) -5A-5A- A)\ t=0 (21) 

where B is the non-Abelian magnetic field constructed from A. The first and last terms 
cancel using Stokes' theorem so we require that A c \ satisfy the Euclidean Yang-Mills 
equations A+V AB+AAB+BAA = with boundary condition A\ t=0 = A. The boundary 
integral Sb[A c i, A] = 0, so we are left with an expression for the small-s dependence of 
#[A S ] 

^[A s ] ^ s ~^I d ^ dttr ( A ci+ B ci) = s a (22) 

which, being a positive power of s, goes to zero with s. This gives a contribution to /(A) 
from the cut along the negative real axis in the vicinity of the origin of 

1 t ds w„ u n sin(ira) A f°° dx At „ , m , 
e A(s-i) s a = v — 1_ e -\ / e ~ Xx x a (23) 



2ixi J s — 1 7i Jo 1 + x 

Since 1/(1 + x) < 1/x for positive x, the last integral is less than X~ a T(a) and we obtain 
the bound 

I 1 f ds e \{s-i) s a | \sin(ira) \ T(a) = 1 

l 27riJ s-1 1 rr\ a e x |T(1 - a) | A°e A ' 1 ' 

which is negligible for sufficently large A. So we conclude that, for large A we can ignore 
this semi-classical contribution to /(A) and reconstruct the vacuum functional as ^[A] pa 
J2 n Tfa+iT ' wnere the ijj n are local functionals of the field computable from a knowledge 
of the vacuum functional evaluated for fields that vary slowly on the scale of the lightest 
glueball mass. 
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